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Abstract 

^ ' A versatile new approach for calculating the depletion potential in a hard 



C^ 



sphere mixture is presented. This is valid for any number of components and 
for arbitrary densities. We describe two different routes to the depletion po- 
tential for the limit in which the density of one component goes to zero. Both 

^5 ' routes can be implemented within density functional theory and simulation. 

Q ' We illustrate the approach by calculating the depletion potential for a big 

I ^, . hard sphere in a fluid of small spheres near a planar hard wall. The density 

functional results are in excellent agreement with simulations. 
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When the separation between two big colloidal particles suspended in a fluid of small 
colloidal particles or non-adsorbing polymers is less than the diameter of the small ones 
exclusion or depletion of the latter occurs leading to anisotropy of the local pressure which 
gives rise to an attractive depletion force between the big particles. Asakura and Oosawa fl]] 
first described this depletion mechanism, suggesting that it would drive phase separation in 
colloid-polymer mixtures. Using excluded volume arguments they calculated the force for 
two hard spheres of radius Rb in a fluid of small hard spheres of radius Rs and showed that 
this is attractive for all separations less than 2Rs and is zero for larger separations. The 
hard sphere model captures the essence of the depletion phenomenon and can be mimicked 
experimentally by suitable choices of colloidal solutions 0,|^ . Depletion forces are of fun- 
damental statistical mechanical interest since they are purely entropic in origin and much 
theoretical @-0 and simulation p,^ effort has been devoted to their determination. Many 
recent experiments, using a wide variety of techniques, have investigated depletion forces 
for colloidal mixtures or for colloid-polymer mixtures 0. Much of the impetus for these 
studies stems from a need to understand how depletion determines phase separation and 
flocculation phenomena and this has stimulated a growing interest in ascertaining quantita- 
tive details of the depletion force. Several of these studies are concerned with a big particle 
near a planar wall but Dinsmore et al. ^ have employed video microscopy to monitor the 
position of a single big colloid immersed in a solution of small colloids inside a vesicle - a 
system resembling hard spheres conflned inside a hard cavity. The experimental conditions 
in Ref. [^ correspond to a small sphere packing fraction rjs = 0.3, which is sufficiently high 
that one expects the depletion force to be substantially different from that given by the 
Asakura-Oosawa result; the latter is valid only in the limit rjs ^ 0. Thus, in order to model 
this and other experimental situations one requires a theory of depletion which is reliable 
at high packings and which can tackle rather general "confining" geometries. The latter 
could represent a planar wall, a wedge or cavity or, indeed, another big particle |jTO|. Here 



we present such a theory for the depletion potential based on a density functional treatment 
(DFT) of a fiuid mixture. Our approach is versatile [|ll|] and avoids the limitations of the 
virial expansion ^ and of the Derjaguin approximation |P,|^ and avoids many of the limi- 
tations inherent in integral-equation theories @,^J^ of the depletion force. It has a distinct 
advantage over a direct implementation of DFT in that it does not require a minimization 
of the free energy in the presence of the big particle. It also suggests an alternative simu- 
lation procedure for calculating depletion potentials. We demonstrate the accuracy of our 
approach by comparing our results for the depletion potential between a big hard sphere 
and a planar hard wall with those of our new simulation and with those of Ref. . 

We consider a multicomponent mixture in which species i {i = 1,2, . . .u) has a chemical 
potential /Xj and is subject to an external potential Vi(r) giving rise to the equilibrium number 
density profiles {pj(r)}. The quantity of interest Wt{ri) = fiui'Ti; {fii},T) — Qti{oo; {fii},T) 
is the difference in grand potential flu between a test particle of species i located at position 
Fj and one located at Tj = cxd, i.e., deep in the bulk, far from the object exerting the external 
potentials. Using the potential distribution theorem in the grand ensemble one can easily 
show fl^ : 



Wt{r,) = -ksT In ( ^^] - l^.(r,) + l^i(oo), (1) 

where T is the temperature. Equation (|l]) is a general result; it is valid for arbitrary inter- 



particle forces. It is important for subsequent application to note that pj(r) is determined 
solely by the external potentials so it has the symmetry dictated by these potentials and 
not the broken symmetry brought about by inserting the test particle. In order to obtain 
the depletion potential W, which pertains to a single (big) particle, i = b, in the presence of 
a fixed object exerting a potential Vf,(r) of finite range, we take the limit fib —>■ — oo so that 
the density of the big particles vanishes, with the other chemical potentials {yUj^fe} fixed. It 
follows that for r outside the range of H(r) 

pW{r) = - lim \n(^^] (2) 

with f3 = [kBT)^^. Although both pfe(r) and pb{oo) vanish in this limit their ratio is non- 
zero. Moreover, the density profiles of the other species {^^^^(r)} reduce to those in the i/ — 1 
component fluid. Thus, for a binary mixture of big and small (s) particles, Ps(r) reduces to 
the profile of a pure fluid of small particles. Equation (^ can be re-stated in a more familiar 
form as p{r)/p{oo) = exp(— /3iy(r)), where p{r) is the probability density of finding the big 
particle at a position r from the fixed object. 

There are two distinct ways of implementing this route to the depletion potential. First 
one can calculate the density profile P6(r) and, hence, Wt{r) in the mixture for decreasing 
concentration of species b. For sufficiently negative values of pb, Wt{r) should approach its 
limiting value W{r). Such a procedure is straightforward to implement in any approximate 
DFT, especially for objects with planar or radial symmetry (e.g., planar walls, spherical 
cavities or a fixed spherical particle). Second one can attempt to proceed immediately to 
the limit of a single big particle, thereby obtaining W{r) directly without evaluation of Pb(r). 
Equation ([^) can be re-expressed as a difference in one-body direct correlation functions |l! 



- (3W,{r,) = cf v.; Ka}]) - cf V; [{P^}]) (3) 

and the depletion potential can be expressed formally as 

- /3Wir) = 4')(r; [{p^^b},Pt = 0]) - 4'V; [{p^^t},p, = 0]). (4) 

In order to implement this result an explicit prescription for c^ is required. Within the 
context of DFT 4^^*; [{Pi}]) = ~/3SJ^ex[{Pi}]/5piiri), where Tex[{Pi}] is the excess (over 



the ideal gas) free energy functional of the inhomogeneous mixture [|r^ . For certain classes 
of (approximate) functionals explicit results can be obtained for c\ which permit the limit 



Pb — ^ — oo to be taken |14]. An important example is Rosenfeld's Fundamental Measure 



functional for hard sphere mixtures |T3 



^ex[{A}]= /"rfV<l>(K(r)}), (5) 

where the Helmholtz free energy density $ is a function of the set of weighted densities 

^"W = E/rfVp,(r')^r(r-r') (6) 



and explicit expressions are available ||T5[ for the weight functions tuf and for $. In this 
theory the free energy of the homogeneous mixture is identical to that from Percus-Yevick 
or scaled-particle theory. Within this approach Eq. (||) reduces to 

pWt{r) = f dV Y^ ^"(r'X(r' - r) (7) 



with i = b and 

i.e., the grand potential difference consists of a sum of convolutions of the functions \&" and 

the weight functions uj" of the big particle. The index a labels 4 scalar plus 2 vector weights 

T^,ITH]. Since the derivatives in Eq. (H) are evaluated at equilibrium, the weighted densities 



must be obtained from Eq. (^, once the density profiles Pj(r) are obtained by minimizing 
the free energy functional. However, for the binary mixture the limit Pb(r) — > implies that 
^^(r) involves only Ps(r) which is given by the solution of the Euler-Lagrange equation for 
the pure small-sphere fluid. The depletion potential W{r) is then given by Eq. (^ with na{r) 
determined in this way. Both schemes have been used to calculate the depletion potential 
for a big hard sphere, immersed in a sea of small spheres, near a planar hard wall. Before 
describing our results we remark upon two limiting cases of the DFT treatment. 

(i) The present functional |jl5|] reduces to the exact low density limit l3J-'ex[{Pi}] = 
~2 Si i / / d^^d^''^' Pi{'^)Pj{'^') fiji'^ ~ ^')y where fij is the Mayer bond between species i and 
j, when all the densities are small and the depletion potential simplifies to 

(3Wir) = -Jd^'ipsir') - p.(oo))/,,(r - r') (9) 

where, in the same limit, Ps(r) = ps(oo) exp(— /5K((r)). For hard spheres /bs(r — i"') = 
—Q{{Rs + Rb) — |r — r'l) where is the Heaviside function and one can show that for a 
hard potential K(r), —l3W{r) is Ps(oo) times the overlap volume of the exclusion sphere 
of the big particle with the hard wall, i.e., one recovers the Asakura-Oosawa result for the 
depletion potential 0]. 

(ii) In the limit where the two species have equal radii, i.e., s = Rs/Rb = 1, it is 
straightforward to show that c^ reduces to the direct correlation function of a pure (s) fluid 
and that the depletion potential is given by j3W{r) = — ln(ps(r)/ps(oo)), which is the exact 
result [ffi. 



Since the Rosenfeld functional is known to yield very accurate results [^ for the density 
profile of a pure fluid near a hard wall and for the radial distribution function, obtained by 
fixing a particle at the origin and calculating the one-body density profile, we have good 
reasons to expect that the depletion potential calculated from DFT for s < 1 should be 
rather accurate for all (fiuid) densities ps(oo). On the other hand for small size ratios 
the accuracy of the mixture DFT is not known; we test this by making comparison with 
simulations. 

In Fig. |I| we show Wt{z) for the binary hard sphere mixture, with s = 0.2, at a planar 
hard wall as calculated using the Rosenfeld functional for a fixed packing fraction rj^ = 



47r_RgPs(oo)/3 = 0.2 of the small spheres and 4 different values of rib = 47r_R^pb(cx))/3. z 
is the distance between the surface of the big sphere and the wall. For rjf, = 10~^, Wt{z) 
is indistinguishable from the depletion potential W{z) calculated using the second, direct, 
method. Moreover, at this concentration Ps{z) is indistinguishable from the the profile of the 
pure fluid. As both functions converge fairly rapidly to their limits (the results for rjb = 10^^ 
differ from those for 10^^ by at most 1.6%) we conclude that this is a practicable method 
of calculating depletion potentials via DFT. We note that for small values of rjb, Wt{z) 
(and Ps{z)) exhibit decaying oscillations whose period is ~ 27?^, whereas for rjb ^ 5 x 10~^ 
additional, longer period oscillations develop. The latter reflect the ordering of the big 
spheres which becomes pronounced at high packings. Figure ^ shows the depletion potential 
W{z) calculated for a size ratio s = 0.2 and rjg = 0.3 alongside the result from Ref. [§. The 
latter was obtained from canonical Monte Carlo simulations of the depletion force, as given 
by an integral of the local density of small spheres in contact with the large sphere 0] . Also 
plotted in Fig. |^ are our own simulation results, based on the formula p{z) oc exp{—j3W{z)). 
As direct measurement of the probability density p{z) of finding the big particle at a distance 
z from the wall yields poor statistics, we measured the probability ratio p{z)/p{z + Az) using 
umbrella sampling |T^ in a grand canonical Monte Carlo simulation of the small spheres, 
i.e., we determined P{W{z + Az) — W{z)), where Az is a step length, and we set W{0) 
equal to the DFT value. Our DFT results are in excellent agreement with those of both 
simulations: the height of the depletion barrier (the maximum value of W minus the contact 
value W^(0)) is given almost exactly and the subsequent extrema are closely reproduced by 
the theory [T^. That two independent sets of simulation results, based on totally different 
routes to the depletion potential, agree so closely is pleasing and attests to the accuracy 
with which these potentials can be calculated. We find a similar high level of agreement 
between DFT and the results of Ref. |^ for rjs = 0.1, 0.2 and the smaller size ratio s = 0.1. 
The oscillatory behavior of W{z) warrants further discussion. From the general theory of 
the asymptotic decay of density profiles in mixtures with short-ranged interaction potentials 
it is known that the profiles of all species near a wall (or a fixed particle) exerting a finite 
ranged external potential exhibit the same characteristic decay [|l8l. Thus, for a binary 
mixture near a hard wall Pi{z) — pj(oo) ~ A^iexp{—aQz) cos{aiz — Qwi), z —^ oo where 
i = b,s. The amplitude A^i and phase 9^j depend on the particular species but a^ and ai 
are common. In the limit pb{oo) -^ these are given by the solution of 1 — ps(oo)css(a) = 0, 
where Css{a) is the Fourier transform of the two-body direct correlation function of the 
homogeneous pure fluid of density ps{oo). ai is the real and Oq the imaginary part of the 



solution a with the smallest imaginary part [|T8|. We have confirmed from our numerical 
results that for r]b = 10~^ and s = 0.2, Pb{z) and Ps{z) exhibit the same value of Oq and ai. 
Moreover, for fixed rjs, both oq and ai are, as predicted, independent of the size ratio s. It 
follows from Eq. (||) that —(3W{z) ~ Au]bexp{—aoz) cos{aiz — Qwb), z ^ oo, and apart from 
the amplitude and phase, the asymptotic decay of the depletion potential is also independent 
of s [0. From a fit to the numerical data for rjg = 0.2 we obtained a^ ^ 2.1/(2/?^) and 
ai ^ 5.2/(2i?s). These values are close to those obtained from solving 1 — ps{oo)css{a) = 
for the pure fluid where, consistent with the Rosenfeld DFT, Css is given by Percus-Yevick 
theory. Moreover, and in keeping with earlier studies of bulk pair correlation functions 
pTSl , we find that the leading-order asymptotic result provides a good description of the 



intermediate as well as the longest-ranged behavior of the depletion potential. We note that 



the period of the oscillations 271 /ai decreases whereas the decay length a^^ increases with 
increasing 77^. 

The success of the comparisons between our DFT results and those of simulation lead 
us to expect that the Rosenfeld functional will provide an accurate description of depletion 
potentials for size ratios down to s = 0.1 and for packing fractions rjg up to 0.3. (For 
smaller size ratios we expect the Percus-Yevick theory of hard sphere mixtures, upon which 
the functional is based, to become inaccurate and there is no reason to expect that any 
approximate mixture functional would reproduce the exact s — > results for the depletion 
potential which are obtained by making the Derjaguin approximation at the outset |^.) 
This implies that our procedure can be profitably employed for some of the more complex 
geometries investigated in experiments 0. We emphasize that our present procedure has 
a crucial advantage over the brute-force application of DFT to the calculation of depletion 
potentials. In the latter one would calculate W{v) either from the total free energy of the 
inhomogeneous fluid or from the local density of the small particles in contact with the big 
one [0,Q. In contrast with the present scheme, both methods require much numerical effort 
to minimize the free energy functional since the original symmetry of the density profile 
Ps(r) is broken by the presence of the big particle. A likely limitation of our procedure is 
that one requires a reliable free energy functional for the (binary) mixture of big and small 



particles. Although accurate functionals are available for hard sphere mixtures 0] this is 
not the case for most other types of mixtures. 

As a final remark we note that it is possible to take simulation data for Ps(r), computed 
in the absence of the big test particle, and insert these into the DFT expression Eq. (^. 
Although such a procedure is not self-consistent, in view of the high accuracy of the Rosen- 
feld functional it does provide an alternative means of calculating depletion potentials for 
complex geometries such as wedges where a direct simulation of the potential or force is 
extremely difficult. 

We are grateful to R. Dickman for providing his simulation results for the depletion 
potential which we plot in Fig. ^. We benefited from conversations with R. Van Roij. This 
research was supported in part by the EPSRC under GR/L89013. 
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FIG. 1. (3Wt{z) calculated for a binary hard sphere mixture with size ratio s = 0.2 and fixed 
small sphere packing fraction r/^ = 0.2 and for 4 different values for r]b. The separation between 
the surface of the big sphere and the hard wall is in units of the small sphere diameter ag = 2Rs . 
For rji) = 10~^, (3Wt{z) is indistinguishable from the depletion potential W{z) obtained from the 
direct method. 
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FIG. 2. The depletion potential VF(z) between a big hard sphere and a planar hard wall for 
a size ratio s = 0.2 and rjg = 0.3 calculated from DFT using the direct method (line) compared 
with results from Ref. [^ (a,0)) shifted as described in Ref. |17|, and with our present simulations 

(b,n). 



